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We study the dynamics of a two-mode Bose-Einstein con- 
densate in the vicinity of a mean-field dynamical instability. 
Convergence to mean-field theory (MFT), with increasing to- 
tal number of particles N, is shown to be logarithmically slow. 
Using a density matrix formalism rather than the conven- 
tional wavefunction methods, we derive an improved set of 
equations of motion for the mean-field plus the fluctuations, 
which goes beyond MFT and provides accurate predictions for 
the leading quantum corrections and the quantum break time. 
We show that the leading quantum corrections appear as de- 
coherence of the reduced single-particle quantum state; we 
also compare this phenomenon to the effects of thermal noise. 
Using the rapid dephasing near an instability, we propose a 
method for the direct measurement of scattering lengths. 



I. INTRODUCTION 

The effective low-energy Hamiltonian for N interacting 
bosons confined in an external potential V ex t , is given in 
second-quantized form as 



H 



2m 



V(r) 
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where V(r) is the external trapping potential, and m is 
the particle mass, g is a coupling constant proportional 
to the s-wave scattering length, and ip,^ are bosonic 
annihilation and destruction operator fields obeying the 

canonical commutation relation ip(r), ^(r') = S(r — 

r'). (This Hamiltonian is an effective low-energy ap- 
proximation, in the sense that short wavelength de- 
grees of freedom have been eliminated: it is applicable 
in the regime of ultracold scattering, where short dis- 
tance modes are only populated virtually, during brief 
two-body collisions.) At very low temperatures, Bose- 
Einstein condensation occurs, so that a large fraction of 
the particles occupy the same single-particle state, char- 
acterized by the single particle wave function ^(r, t). In 
this regime one can formulate a perturbative expansion 
in the small quantity A -1 / 2 , where N is the number of 
particles in the condensate, whose result at leading order 
is the Gross-Pitaevskii nonlinear Schrodinger equation 
(GPE) governing the condensate wave function: 
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The Gross-Pitaevskii mean field theory (MFT) pro- 
vides a classical field equation for nonlinear matter waves, 
which is generally considered as 'the classical limit' of 
the Heisenberg equation of motion for the field opera- 
tor ip (which is of precisely the same form). We can 
make precise the sense in which it is a classical limit, 
by reformulating the system governed by (^) in the path 
integral representation. We will not actually use this 
formulation in this paper; we merely note the GPE is the 
saddlepoint equation that appears in a steepest descents 
approximation to the path integral. This is precisely 
the standard semi-classical approximation, with the ex- 
ception that l/N is playing the role usually played by 
%. Hence despite the resemblance of the GPE to a 
Schrodinger equation, complete with finite h, we can in- 
deed identify MFT as the classical limit, in essentially 
the same sense as in the case Ti — > 0, of the quantum field 
theory. Because N in current trapped dilute alkali gas 
BEC experiments is characteristically large (typically of 
the order I0 5 — 10 s atoms), qualitatively significant quan- 
tum corrections to MFT are hard to observe, and the GP 
theory is highly successful in predicting experimental re- 
sults. 

The entire field of quantum chaos is founded upon one 
property of the classical limit, however, which is that 
convergence to classicality as h — ► is logarithmically 
slow if classical trajectories diverge exponentially. This 
implies that we must expect strong quantum corrections 
to MFT in the vicinity of a dynamically unstable fixed 
point. In particular, the quantum evolution will de- 
part significantly from the classical approximation af- 
ter a logarithmic 'quantum break time', which will be 
"log A in our case, as it is ~log(l/7i) in the standard 
case . In our case, the nature of this departure is 
that after the quantum break time, a condensate will 
become significantly depleted, as exponential production 
of quasi-particles transfers particles to orthogonal modes 
PJ. Depletion of the condensate means, by definition, 
that the single particle reduced density matrix (SPDM) 
becomes quantum mechanically less pure. Hence for a 
condensate, just as the classical limit of the quantum 
field theory resembles the quantum mechanics of a single 
particle, so quantum corrections at the field theory level 
appear as quantum decoherence in the single-particle pic- 
ture. Since decoherence is most often considered as en- 
forcing classicality, there is something like irony in this 
situation. And it suggests that studying the corrections 
to MFT for Bose-Einstein condensates may give us some 
new insights into decoherence; and that some aspects of 
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decoherence may be useful in understanding condensates 
beyond MFT. This is the motivation for the work we 
now report. 

In this paper we provide the details of a previously pub- 
lished study |^| of the correspondence between mean-field 
and exact quantum dynamics of a two- mode BEC. The 
model system contains an isolated dynamical instabil- 
ity for certain regions of parameter space. We show that 
quantum corrections in the vicinity of this unstable state, 
do indeed become significant on a short log(iV) time 
scale, whereas quantum effects in other regions of phase 
space remain small 1/yN corrections. We present a sim- 
ple theory that goes beyond MFT and provides accurate 
predictions of the leading quantum corrections, by taking 
one further step in the so called Bogoliubov-Born-Green- 
Kirkwood-Yvon (BBGKY) hierarchy. In accordance with 
our view of quantum corrections as decoherence, we use 
a density-matrix Bloch picture to depict the dephasing 
process. The density-matrix formalism has the addi- 
tional advantage of allowing for initial conditions that 
are not covered by the Hartree-Fock-Bogoliubov Gaus- 
sian ansatz, and which better correspond to the physical 
state of the system. 
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FIG. 1. Two- mode Bose-Einstein condensates: (a) a con- 
densate in a double-well potential (b) a spinor condensate. 

In section II we briefly review the model system and 
its experimental realizations. In section III we derive the 
mean-field equations of motion in the Bloch representa- 
tion, and illustrate the main features of the produced dy- 
namics for various parameter sets. Quantum corrections 
to the two-mode MFT are studied section IV, as well 
as an improved theory that predicts the leading correc- 
tions. In section V we consider the effect of thermal noise, 
and show an analogy between the quantum dephasing of 
the reduced single particle density operator and thermal 
decoherence. In section VI we present a potential ap- 



plication of the rapid decoherence near the dynamical 
instability of the two mode model, for the measurement 
of s-wave scattering-lengths. Discussion and conclusions 
are presented in section VII. 



II. THE TWO-MODE CONDENSATE 

We consider a BEC in which particles can only ef- 
fectively populate either one of two second-quantized 
modes. Two possible experimental realizations of this 
model are illustrated in Fig. 1. The first (Fig la) is 
a condensate confined in a double-well trap |3|-|£]] which 
may be formed by splitting a harmonic trap with a far 
off resonance intense laser sheet || . In this case single- 
particle tunneling provides a linear coupling between the 
local mode solutions of the individual wells, which can 
in principle be tuned over a wide range of strengths by 
adjusting the laser sheet intensity. The two-mode regime 
is reached when the self-interaction energy gn is small 
compared to the spacing between the trap modes hcotrap- 
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where I is the characteristic trap size. Thus the two-mode 
condition is, 



I > N\a\ 



(4) 



The two-mode condition (^) may be met by double-well 
traps with characteristic frequencies of the order of 100 
Hz, containing several hundred particles. When con- 
structed, larger traps will maintain the two-mode limit 
at higher N. 

The second experimental realization of a two-mode 
BEC is the effectively two-component spinor condensate 
0,0 depicted in Fig. lb. In this case the linear cou- 
pling between the modes is provided by a near resonant 
radiation field |I^|l3|]. If collisions do not change spin 
states, the nonlinear interactions between the particles 
depend on three scattering lengths a^ . In realisations 
of spinor condensates it is easy to ensure an = 022 by 
symmetry, in which case the nonlinear interaction term 
becomes H int = Jd 3 rH int (r) for 
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We can therefore define two healing lengths £± = 
\/ yj p{a\\ ± 012), where p is the total density, charac- 
terizing the effect on the spatial state of the condensate 
of the two nonlinear terms. The two-mode regime, in 
which the spatial state is fixed and essentially indepen- 
dent of the internal state, is reached when £_ becomes 
larger than the sample size (its largest dimension) . Since 
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for available alkali gases all differ only by a few per- 
cent, << £+, and hence the two-mode regime can be 
reached with N < 10 4 atoms in weak, nearly spherical 
traps (cotrap < 100 Hz). Less isotropic traps obviously 
reach the two-mode regime only at smaller N. To ex- 
tend the internal state two mode regime to larger N we 
must make the trap weaker; for fixed N, £_ scales with 
sample size L as L 3 / 2 . Hence to ensure < 1 for 

fixed N requires a sufficiently large (weak) trap. For 
Rb and Na experiments, whose lifetimes are limited by 
three-body collisions, the slowing down of the two-mode 
dynamics at reduced total density should be more than 
compensated for by the extended condensate lifespan. 

In both realizations, the many-body Hamiltonian re- 
duces in the two-mode limit (and in the spinor realization 
also in the rotating-wave approximation) to the form, 



H(t) = {a[ai + a 2 a 2 ) - — [a[ 
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where E\ and E% are the two condensate mode ener- 
gies, Q is the coupling strength between the modes, g is 
the two-body interaction strength, and ai,a[, ci2,a\ are 
particle annihilation and creation operators for the two 
modes. The total number operator N = a\a\ + a\a,2 
commuted with H and may be replaced with the c- 
number N. Writing the self-interaction operators as 
(a\) 2 a\+{a\) 2 a\ = [N 2 +{d\ai—a\a.2) 2 }/2 and discarding 
c-number terms, we obtain the two-mode Hamiltonian 
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we obtain 



h = -mt x + —L 2 
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The Heisenberg equations of motion for the three angular 
momentum operators of Eq. (pf) read 
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— I/a; — —~[L X , H] — —~(LyL Z + L z Ly) . 

^L y = -^[L y ,H]=+QL z + ^(L x L z + L z L x ) , (10) 
d - i « 

-L z = --[L z ,H\ = -m y . 

Thus the expectation values of the first-order opera- 
tors Li depend not only on themselves, but also on the 
second-order moments (LiLj). Similarly, the time evo- 
lution of the second-order moments depends on third- 
order moments; and so on. Consequently, we obtain 
the BBGKY hierarchy of equations of motion for the 
expectation- values, 



dt 
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(L i L j )=fl(L i ,L j ,),(L i ,L j ,L k ,)) , (11) 
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We will take g and ui to be positive, since the relative 
phase between the two modes may be re-defined arbi- 
trarily, and since without dissipation the overall sign of 
H is insignificant. 



III. TWO-MODE MEAN-FIELD THEORY IN THE 
BLOCH REPRESENTATION 



where i, j,k, . . . , i , j , k' , I', . . . — x, y, z. In order to ob- 
tain a closed set of equations of motion, the hierarchy of 
Eq. (11) must be truncated at some stage by approxi- 
mating the iV-th order expectation value in terms of all 
lower-order moments. 

The lowest-order truncation of Eq. ( pi] ) is obtained 
by approximating the second-order expectation values 
(LiLj) as products of the first-order moments (Li) and 
(L<): 



The conventional wavefunction formalisms consider 
the evolution of dj and its expectation value in a 
symmetry-breaking ansatz (where the symmetry being 
broken is that associated with conservation of N). In- 
stead, we will examine the evolution of the directly ob- 
servable quantities ajoj, whose expectation values de- 
fine the reduced single particle density matrix (SPDM) 
Rij = (d\a,j)/N . Writing the Hamiltonian of Eq. (0) in 
terms of the SU(2) generators, 
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The equations of motion for the single-particle Bloch vec- 
tor 



Then read 
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where k — gN/2. Equations ([[J]) describe rotations of 
the Bloch vector s, and so the norm |s| is conserved 
in MFT. Consequently, for a pure SPDM, Eq. (|J) are 
completely equivalent to the two- mode Gross-Pitaevskii 
equation ra], 
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where a± and 02 are the c-number coefficients replacing 
the creation and annihilation operators of Eq. (ph 
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FIG. 2. Mean-field trajectories at (a)n = 0, (d)k = 1.02O, 
(c)k = 20., and (d)n = 20O. 

In Fig. 2 we plot mean-field trajectories at four differ- 
ent n/Q ratios. The nonlinear Bloch equations (|l4| ) de- 
pict a competition between linear Rabi oscillations in the 
Sj,s z -plane and nonlinear oscillations in the s^s^-plane. 
For a noninteracting condensate (Fig. 2a) the trajecto- 
ries on the Bloch sphere depict harmonic Rabi oscillations 
about the s x axis. As n increases the oscillations become 
increasingly anharmonic. As long as k < Q the nonlincar- 
ities may be treated as perturbation. However, above the 
critical value k = uj (Fig. 2b), there are certain regions 
in phase-space which are dominated by the nonlinear 
term. The stationary point s = (—1,0,0), correspond- 
ing to the Josephson 7r-state (equal populations and a 7r 
phase-difference), becomes dynamically unstable and the 
two trajectories passing asymptotically close to it form 
a "figure-eight" . The region outside these limiting tra- 
jectories is dominated by the linear oscillations whereas 
inside the nonlinear term prevails. Starting at the criti- 
cal value of k = 2f2 (Fig. 2c) population prepared in one 
of the modes remains trapped in the half-sphere it origi- 
nated from, conducting oscillations with a non- vanishing 
time averaged population imbalance (s z }t ^ 0. This phe- 
nomenon was termed "macroscopic self-trapping" H . Fi- 



nally, when k 3> Q (Fig. 2d) the nonlinearity dominates 
the entire Bloch sphere, except for a narrow band about 
the s z = plain. 



IV. QUANTUM CORRECTIONS AND 
BOGOLIUBOV BACKREACTION 

In the vicinity of the dynamically unstable point, we 
expect MFT to break down on a time scale only log- 
arithmic in TV. In order to verify this prediction, we 
solve the full TV-body problem exactly, by fixing the total 
number of particles TV, thereby restricting the available 
phase-space to Fock states of the type | n, N — n ) with 
n particles in one mode and TV — n particles in the other 
mode, n ranging from to TV. Thus we obtain an TV + 1 
dimensional representation for the Hamiltonian (|^) and 
the TV-body density operator p: 



H m ,n = (m, TV — m| H \ n, TV — n) 
p m ,n = (m,N -m\p\n,N -n) , 



(16) 
(17) 



for m,n — 0,1,...,TV. The exact quantum solution is 
obtained numerically by propagating p according to the 
Liouville von-Neumann equation 



ihp = [H, p] 



(18) 



Using the Hamiltonian of Eq. (|9|) to evaluate the matrix 
elements of Eq. ([l6]) and substituting into Eq. (filj|), 
we obtain dynamical equations for the N-body density 
matrix: 

n 



ihp m ,n = -— \/ m(N -m + l)p m -i, n 



+ \/(m + 1) (TV - m)pm+ 1 , n 
- \Jn(N -n + l)p m , n _i 



(19) 



- yin + 1)(TV - n)p m . n+ i 



+ 7 [m 2 — (TV — mf - n 2 + (TV - n) 2 ] p n 



Equations ( p^ ) are solved numerically, using a Runge- 
Kutta algorithm. In fig. 3 we plot exact quantum trajec- 
tories starting with all particles in one mode, for increas- 
ingly large TV (k being fixed) versus the corresponding 
mean-field trajectory. While MFT assumes a persistently 
pure single particle state, quantum corrections to MFT 
appear in the single-particle picture, as decoherence of 
the SPDM. When the mean-field trajectory stays away 
from the instability (Fig. 3a) the quantum trajectories 
indeed enter the interior of the unit Bloch sphere at a rate 
that vanishes as l/\/TV. However, when the mean-field 
trajectory includes the unstable state (Fig. 3b), we ob- 
serve a sharp break of the quantum dynamics from the 
mean-field trajectory at a time that only grows slowly 
with TV. 
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In accordance with our picture of quantum corrections 
as decoherence, and in order to obtain a more quantita- 
tive view of the entanglement-induced dephasing process, 
we plot the von Neumann entropy 



S = Tr(R\nR) = --In 



(l + |g|)(i+l*l)(l_|g|)(i- 
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of the exact reduced single-particle density operator, as 
a function of the rescaled time ftt for the same initial 
conditions as in Fig. 3. The results are shown in Fig. 
4. Since the entropy of mean-field trajectories is identi- 
cally zero, S may serve as a measure of the deviation from 
MFT. When the mean-field trajectory is stable (Fig. 4a), 
the single particle entropy grows at a steady rate which 
vanishes as N is increased. The variations in the en- 
tropy growth curve are a function of the distance from 
the instability. Near the instability (Fig. 4b) quantum 
corrections grow rapidly at a rate which is independent of 
N, and the time at which this divergence takes place (the 
quantum break time) evidently grows only as log(iV) . 



FIG. 3. Exact quantum trajectories starting with all par- 
ticles in one mode, with N = 50 ( ), 100 ( ), 200 (■ ■ ■) 

and 400 (— • — ) particles vs. the corresponding mean-field 
trajectory ( ) for (a) K = fi and (b) k — 2Q. 

Since MFT can thus easily fail near dynamical insta- 
bilities, it is highly desirable to obtain an improved the- 
ory in which Bloch-space trajectories would be allowed to 
penetrate into the unit sphere without having to simulate 
the entire iV-body dynamics. In fact, such an improved 
non-unitary theory is easily derived using the next level 
of the BBGKY hierarchy (|i"T|). This hierarchy trunca- 
tion approach is in fact a systematic perturbative ap- 
proximation; but it is state-dependent. That is, it pro- 
vides a perturbative approximation, not to the general 
evolution, but to the evolution of a special class of ini- 
tial states, within which the perturbative parameter is 
small. In the case of ultracold bosons, the phenomenon of 
Bose-Einstein condensation ensures that there is a com- 
monly realisable class of states in which the system is a 
mildly fragmented condensate. In our two mode model, 
this means that the two eigenvalues of R are / and 1- f 
for /< 1; and so from such initial states we can approx- 
imate the evolution perturbatively using / as our small 
parameter. 
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To zeroth order in /, R is by definition a pure state, 
and hence we have the MFT evolution on the surface 
of the Bloch sphere. Going to next order in / can be 
achieved by truncating the BBGKY hierarchy at one or- 
der higher. We take Li = Li + 8Li, where the c- number 
Li is O(N) and all the matrix elements of 5 Li remain 
smaller than 0{Ny/J) throughout the evolution of the 
system. The second order moments, 
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will then be of order /. Writing the Heisenberg equa- 
tions of motion for the first- and second-order operators 
Li,LiLj, taking their expectation values and truncating 
Eq. (|ll|) by approximating 



{LiL 3 L k ) » (LiL 3 )(L k ) + (L % }(L 3 L k ) + (LiL k )(Lj) 



2(L l )(L j )(L k ) 
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instead of the mean- field approximation (|12|), we ob- 
tain the following set of nine equations for the first- and 
second-order moments: 
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Equations ( J23| ) will be referred to as the "Bogoliubov 
backreaction equations" (BBR), because they demon- 
strate how the mean-field Bloch vector s drives the fluc- 
tuations Aij - which is the physics described by the 
Bogoliubov theory of linearized quantum corrections to 
MFT; but they also make the Bloch vector subject in 
turn to backreaction from the fluctuations, via the cou- 
pling terms —K,A yz /2 and kA xz /2. This back-reaction 
has the effect of breaking the unitarity of the mean-field 
dynamics. Consequently, the BBR trajectories are no 
longer confined to the surface of the Bloch sphere, but 
penetrate to the interior (representing mixed-state Rij, 
with two non-zero eigenvalues). (Obviously, if the trajec- 
tories penetrate the sphere too deeply, so that the smaller 
eigenvalue / ceases to be small, the entire approach of 
perturbing in / will break down.) 

In order to demonstrate how the BBR equations (|22 
improve on MFT, we compare trajectories obtained by 
these two formalisms to the exact 50-particle trajectories 



of Fig. 3. Both the k = fl stable mean-field trajectory 
and the k = 2fl unstable mean-field trajectory cases are 
plotted in Fig. 5a and Fig. 5b, respectively. The initial 
conditions for the BBR equations are determined by the 
initial state \N,0) to be 



Sz = -1, 

A XX = Ayy = 2/N, 

S x — Sy — ^xy — ^xz — ^-yz — ^-zz 



(24) 
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FIG. 5. Mean field (■■■), Bogoliubov back-reaction ( ) 

and exact 50 particles ( ) trajectories starting with all par- 
ticles in one mode, at (a) n = Q and (b) n — 2£l. 

The approximation of Eq. ( |22| ) ignores terms smaller 
than C(/ 3 / 2 ). It is therefore better than the mean-field 
approximation ( |l2| ) by a factor of f 1 ^ 2 . Consequently, as 
is clearly evident from Fig. 5a, the BBR equations (22) 
are far more successful than the mean- field equations (14) 
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in tracing the full quantum dynamics. However, for any 
realistic number of particles, the improvement is hardly 
necessary, as MFT would be accurate for very long times. 
On the other hand, when the mean-field trajectory ap- 
proaches the instability (Fig. 5b), the BBR theory pro- 
vides an accurate prediction of the leading quantum cor- 
rections. Of course, since the BBR equations account for 
only six moments unobserved by MFT, the period of the 
still quasiperiodic BBR motion is shorter than that of 
the exact evolution and the BBR trajectory eventually 
deviates from the quantum trajectory. Nevertheless, the 
BBR formalism provides a simple and efficient method to 
predict the quantum break time in large (> 1000 parti- 
cles) condensates, for which full iV-body simulations are 
restricted by available c omp utation power. 

The BBR equations (G3) are in fact identical to the 
equations of motion one would obtain, for the same 
quantities, using the Hartree-Fock-Bogoliubov Gaussian 
ansatz, in which second-order moments A, 3 - are initially 
factorized as = Sidj — x, y, z). Using this ansatz, 
the factorization persists and the time evolution of 6 x ,d y , 
and 5 Z is equivalent to that of perturbations of the mean- 
field equations (|14|): 



S x = -K(S z 5y + SyS z ) , 

Sy = loS z + k(s z 6 x + s x S z ) , 
<L = —<jj6„ . 



(25) 



Thus our equations for A<y are in a sense equivalent to the 
usual Bogoliubov equations. The quantitative advantage 
of our approach therefore lies entirely in the wider range 
of initial conditions that it admits, which may more accu- 
rately represent the exact initial conditions. For instance, 
a Gaussian approximation will have A xx = 0(1) in the 
ground state, where in fact A xx = C(iV _1 ). This leads to 
an error of order TV -1 / 2 in the Josephson frequency com- 
puted by linearizing ( p3| ) around the ground state, even 
though the Gaussian backreaction result should naively 
be accurate at this order. Our SPDM approach does not 
have this flaw, which is presumably the two-mode version 
of the Hartree-Fock-Bogoliubov spectral gap [B . 



V. DEPHASING DUE TO THERMAL NOISE 

Decoherence is generally considered as suppressing 
quantum effects fll5[ . Ironically, in our case the lead- 
ing quantum corrections to the effectively classical MFT, 
are themselves decoherence of the single-particle state of 
the condensate. Therefore, it is interesting to study the 
effect of a realistic decoherence process, originating in the 
coupling to a bath of unobserved degrees of freedom, on 
the interparticle entanglement process, described in the 
previous section. 

The main source of decoherence in BEC's is the ther- 
mal cloud of particles surrounding the condensate. Ther- 
mal particles scattering off the condensate mean field will 



for example, cause phase diffusion [[16| at a rate T pro- 
portional to the thermal cloud temperature. For internal 
states not entangled with the condensate spatial state, T 
may be as low as 10~ 5 Hz under the coldest experimen- 
tal conditions, whereas for a double well the rate may 
reach 10 _1 Hz. Further sources of decoherence may be 
described phcnomenologically with a larger T. 

We account effect of thermal noise on the two-mode 
dynamics by using the quantum kinetic master equation 
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(26) 



Once again, we solve for s(t) using cither one of three 
methods: 

(a) MFT - The decoherence term in Eq. ( pBf ) introduces 
an exact T2 = 1/r transversal relaxation term into the 
mean-field equations of motion: 



K 1 S Z S- L 



UJS Z + KS Z S 

— UJSy . 



Fs x , 



(27) 



(b) BBR - Evolving the first- and second-order oper- 
ators according to Eq. (p6f), taking their expectation 
values and truncating the hierarchy at the next level, we 
obtain the modified BBR equations, 
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(c) Exact quantum solution - obtained by numerically 
propagating the full A^-particlc density matrix under Eq. 

In Fig. 6 we compare the Von-Neumann entropy of 
the exact TV-body density operator as a function of time 
for exponentially increasing N , to the mean-field entropy. 
Due to the thermal noise, mean-field trajectories are no 
longer confined to the zero-entropy sphere. However, 
whereas the quantum break time in the absence of ther- 
mal noise has grown as log (AT) (see Fig. 4b), it is clear 
from Fig. 6 that in the presence of this dephasing mech- 
anism it saturates to a finite value. Thus, while we may 



7 



have naively expected decoherence to reduce quantum 
corrections and thereby improve MFT, in fact the ad- 
dition of thermal dephasing has significantly damaged 
classical-quantum correspondence. 



0.7 




reseated time 



FIG. 6. Growth of the von Neumann entropy 5* of the quan- 
tum reduced single-particle density operator in the presence 

of thermal noise ( Y = 10~ 2 ft), at K = 20, for N=10 ( ), 

20 (■ ■ ■), 40 ( ), 80 ( ), 160 (- . - ), 320 ( , 

and 640 (. . .) particles. Bold solid curve corresponds to the 
mean field entropy. Initial conditions are the same as in Fig. 
3. 

In Fig. 7, we summarize the results of numerous dy- 
namical calculations conducted for various values of the 
particle number N and of the thermal noise T, by plotting 
the time at which the entropy reaches a given value. The 
curves are obtained using the modified BBR equations 
( p8| ) whereas the circles and squares depict exact quan- 
tum results (limited by computation power to N ~ 10 3 
particles) for two limiting values of T. The BBR equa- 
tions provide accurate predictions of the initial decoher- 
ence rate and the quantum break time even within this 
limited range of N (and the agreement between the exact 
quantum results and the BBR predictions would become 
still better for higher numbers of particles). Once more, 
we observe the logarithmic growth of the quantum break 
time with N in the zero temperature (r = 0) limit. How- 
ever, when the temperature is finite, there is a saturation 
of the quantum break time to values which are well be- 
low the mean-field thermal dephasing times, in agreement 
with Fig. 6. 

Instead of observing the quantum break time as a func- 
tion of the number of particles for a given degree of ther- 
mal noise, we can monitor the thermal decoherence time 
as a function of temperature, for any given number of 
particles. Viewing Fig. 7 this way, it is evident that 
in the mean-field limit (l/y/~N — ► 0) the purely thermal 
dephasing time also grows only logarithmically with the 
temperature, Comparing this result to the log(A) growth 



of the quantum break time in the zero-temperature limit, 
we can see that thermal noise and quantum noise have es- 
sentially similar effects on the system. And Figs. 6 and 
7 together are in complete agreement with the prediction 
that the entropy of a dynamically unstable quantum sys- 
tem coupled to a reservoir |]l7j] , or of a stable system 
coupled to a dynamically unstable reservoir, will grow 
linearly with time, at a rate independent of the system- 
reservoir coupling, after an onset time proportional to the 
logarithm of the coupling pqjig|] . Thus, one can really 
consider the Bogoliubov fluctuations as a reservoir p0[ , 
coupled to the mean field with a strength proportional to 
1 /N. The N «-> 1/T analogy is even further extended by 
the saturation for any finite N, of the thermal dephas- 
ing time at low T, in the same way that the quantum 
break time for a finite T saturates at high N. Due to 
this quantum saturation, quantum corrections can be ex- 
perimentally distinguished from ordinary thermal effects 
which do not saturate the dephasing rate at low temper- 
ature. 









1-0 




y< r=l (T 4 a> 




r=i(T 3 co 
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FIG. 7. Time at which S reaches 0.2 as a function of 
the particle number N, according to the BBR equations (p3|), 
modified to include thermal phase-diffusion. Four different 

values of V are shown: T = ( ), V = 10~ 4 S1 (■■■), 

T = 10~ 3 fi ( ), and V = 10~ 2 £1 ( ). Exact quan- 
tum results are presented for T = (circles) and F = 10 _2 f2 
(squares). Initial conditions, re and Q are the same as in Fig. 
6. 



VI. SCATTERING-LENGTH MEASUREMENTS 

After indicating how condensate decoherence at dy- 
namical instabilities can connect principles established 
in different areas of physics, we briefly note that it can 
also have practical applications. Rapid decoherence in 
the vicinity of the unstable 7r-state of the two-mode con- 
densate may serve for the direct measurement of scatter- 
ing lengths. As demonstrated in Fig. 8, the mean-field 
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trajectory of a condensate which is prepared initially in 
one of the modes, would only pass through the rapidly 
dephasing unstable point when k = 20. Thus, the self- 
interactions energy k can be determined by measuring 
the entropy at a fixed time as a function of the coupling 
frequency f2, resulting in a sharp line about = /c/2, as 
depicted in Fig. 9. 

Experimentally, the single-particle entropy is measur- 
able, in the internal state realization of our model, by 
applying a fast Rabi pulse and measuring the amplitude 
of the ensuing Rabi oscillations, which is proportional to 
the Bloch vector length |s|. (Successive measurements 
with Rabi rotations about different axes, i.e. by two res- 
onant pulses differing by a phase of tt/2, will control for 
the dependence on the angle of s) . In a double well real- 
ization, one could determine the single-particle entropy 
by lowering the potential barrier, at a moment when the 
populations on each side were predicted to be equal, to 
let the two parts of the condensate interfere. The fringe 
visibility would then be proportional to |s| M. 




K — 4Q ( ). The dot at s = ( — 1, 0, 0) marks the dynamical 

instability. 

VII. CONCLUSIONS 

To conclude, we have shown that significant quantum 
corrections to the Gross-Pitaevskii MFT, in the vicin- 
ity of its dynamical instabilities, can be measured in a 
two-mode BEC under currently achievable experimental 
conditions. We have derived a simple theory that accu- 
rately predicts the leading quantum corrections and the 



quantum break time. By applying to condensate physics 
some insights from studies of decoherence, we have found 
evidence that MFT dynamical instabilities cause linear 
growth of the single-particle entropy at a rate indepen- 
dent of N. And from condensate physics we have learned 
something about decoherence: we have identified a form 
of decoherence which degrades quantum-classical corre- 
spondence, instead of improving it. 

0.6 1 1 1 



0.4 



0.2- 




0)(Hz) 

FIG. 9. von Neumann entropy S after 20 ms of propagation 
according to the BBR equations ( pq ) with k = 1 KHz and 
r = 1CF 4 Hz, starting with the entire condensate of N = 10 5 
particles in one mode, as a function of the coupling frequency 

n. 

Our picture of quantum backreaction in BECs as de- 
coherence suggests new lines of investigation for both 
experiment and theory: measurements of single-particle 
entropy in condensates, descriptions of condensates with 
mixed single particle states (instead of the usual macro- 
scopic wave functions), and general questions of deco- 
herence under nonlinear evolution. Exploring these pos- 
sibilities, beyond the two-mode model considered here, 
provides many goals for further research. 
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